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RIGID  ROTOR  EQUILIBRIUM 
OF  A STRONG  ION  LAYER 

I.  INTRODUCTION 

The  conceptual  simplicity  of  rigid  rotor  equilibria  make  them  very 
attractive  models  for  the  theoretical  studies  of  charged  particles  lay- 
ers. Such  rigid  rotor  equilibria  within  the  framework  of  the  Maxwell- 
Vlasov  equations  (DAVIDSON  et  al.,  I969,  1570  and  I97I*  ; ANDERSON  et  al., 
1971  ; MARX,  I968)  have  been  applied  in  recent  years  to  magnetically 
confined  nonneutral  plasmas  and  to  electron  layers  associated  with  the 
Astron  reactor  concept. 

In  the  nonneutral  plasmas  case,  the  model  (DAVIDSON  et  al.,1969, 

1970  and  197^)  was  based  on  the  nonrelativistic  limit  and  the  axial 
diamagnetic  field  as  well  as  the  effect  of  the  conducting  wall  sur- 
rounding the  gyrating  electrons  has  been  neglected.  By  contrast,  the 
diamagnetic  field  was  included  in  the  theoretical  model  studied  in  con- 
nection with  the  E-layer  (ANDERSON  et  al. , I97I  ; MARX,  I968),  but  the 
assumed  distribution  function  did  not  allow  radial  confinement  of  the 
E- layer. 

In  this  paper,  we  make  use  of  the  steady  state  (9/3t  - o) , Vlasov- 
Maxwell  equations  to  study  the  rigid  rotor  equilibrium  properties  of  a 
space  charge  neutral  proton  layer  (p- layer),  that  is  confined  radially 
by  an  externally  applied  magnetic  field.  The  main  difference  between 
the  present  work  and  that  previously  carried  out  in  connection  to  E- 
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layers  is  in  the  choice  of  the  distribution  function.  The  selected  dis- 
tribution function  (CHU  and  KAPEIANAKOS,  1975)  in  the  present  work  not 
only  allows  radial  confinement  of  the  gyrating  protons,  but  also  is 
simple  enough  that  explicit,  analytic  expressions  can  be  obtained  for 
the  particle  density,  current  density,  magnetic  field  and  the  radii  of 
the  p- layer,  without  any  a priori  assumptions  about  the  radial  dimen- 
sions of  the  layer. 

A brief  description  of  the  equilibrium  configuration  and  an  outline 
of  the  general  assumptions  are  given  in  Sec.  II.  The  equilibrium  prop- 
erties are  calculated  for  the  specific  choice  of  the  proton  distribu- 
tion function  f°  (H-cdPq)  = (mn/2rr)  6 (H-ajpg-kj.),  where  m,  n,  cu  and  k± 
are  constants,  H is  the  total  energy  and  Pg  is  the  canonical  angular 
momentum.  In  Sec.  Ill,  we  review  briefly  results  related  to  the  sta- 
bility of  the  p- layer  and  in  Sec.  IV  we  discuss  the  relevance  of  the 
present  equilibrium  to  pertinent  experiments. 


II.  EQUILIBRIUM  CONFIGURATION  AND  ASSUMPTIONS 


The  equilibrium  configuration  is  shewn  schematically  in  Fig.  1.  It 
consists  of  a space  charge  neutralized,  infinitely  long,  proton  layer 
confined  radially  by  an  externally  applied  magnetic  field  Bq.  We  in- 

| . 

troduce  a cylindrical  polar  coordinate  system  (r,  6,  z) , with  the  z 
axis  along  the  axis  of  symmetry.  The  p- layer  rotates  about  the  axis  of 
symmetry  with  a mean  azimuthal  velocity  Vg(r)  in  the  negative  S-direc- 
tion. 

To  make  the  theoretical  analysis  tractable,  the  following  simpli- 
fying assumption  are  made  in  describing  the  p- layer  equilibrium  by  the 
s teady- state  0 Vlasov- Maxwell  equations:  t •' 
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(a)  Equilibrium  properties  are  azimuthally  symmetric  and  indepen- 
dent of  2- coordinate. 

(b)  The  equilibrium  electric  field  of  the  p- layer  is  neutralized  ly 
an  electron  background  and  the  current  carried  by  the  background 
electrons  is  equal  to  zero.  Thus,  the  self  magnetic  field  is 
generated  entirely  by  the  p- layer. 

(c)  The  magnetic  flux  enclosed  by  the  conducting  cylinder  is  con- 
served for  the  entire  life  of  the  layer,  and 

(d)  The  equilibrium  distribution  function  describing  the  p- layer 
does  not  depend  upon  the  axial  velocity  of  protons,  i.e. , there  is 
not  motion  along  the  z-axis. 

It  is  well  known,  that  any  distribution  function  that  is  a func- 
tion only  of  the  single  particle  constants  of  the  motion  satisfies  the 
steady-state  Vlasov  equation.  For  a time  independent,  azimuthally 


mrv, 


are  the  two  constants  of  the  motion 


proton  mass,  q = e is  the  proton  charge,  c the  speed  of  light  in 
vacuum,  v^,  Vq  and  v z are  the  components  of  proton  velocity  and  Ag(r) 
is  the  equilibrium  magnetic  vector  potential.  Note  that  Ag(r)  = Ag‘  (r) 

Cfll  f pv  f* 

+ Ag  (r),  where  Ag'*  (r)  describes  the  externally  applied  magnetic 
field  Bo  and  AgC^(r)  describes  the  axial  self  magnetic  field  B*(r)  of 
the  layer. 


Experimentally,  it  is  desirable  to  form  p-layers  without  shear  in 
their  mean  azimuthal  velocity  of  rotation  and  without  all  the  particles 
composing  the  layer  having  the  same  energy  and  the  same  canonical 
angular  momentum.  Such  equilibria  are  in  general  more  stable  than 
those  having  shear  in  their  mean  azimuthal  velocity  and  no  energy  or 
momentum  spread.  For  these  reasons,  we  are  considering  an  equilibrium 
distribution  function  than  depends  on  H and  Pg  through  the  linear  com- 
bination, H - (jdP g , where  iu  is  a constant.  Specifically,  the  distribu- 
tion function  is  assumed  to  be  of  the  form 

/°(H-cuPq)  = (mn/2ri)  6 (H-cuPg-ki),  (5) 

where  m,  n and  kj.  are  constants.  Using  Eq.  (1)  with  v = o and  Eq.  (2) 

z 

the  argument  of  the  delta  function  in  Eq.  (3)  may  be  expressed  as 
H - u)Pg  - « (m/2)  [(vg-ru))2  + v^]  + U(r),  (4) 

where 

U(r)  « - (qr/c)  A^(r)cu  - mr 2uj2/2  - kx  . (5) 

The  equilibrium  density  profile  corresponding  to  the  distribution 
of  Eq.  (3)  can  be  computed  from 

n°(r)  =■  JY?(H-uiPg)d2v.  (6) 

Substituting  Eqs.  (3),  (4)  and  (5)  into  Eq.  (6),  it  is  straightforward 
to  show  that  the  proton  density  profile  is 

o,  o ^ r < ai  , 

n°(r)  = n,  ax  * r s a2,  (7) 

o,  a.z  < r < b , 

where  ai  and  a2  are  the  roots  of  function  U(r)  defined  in  Eq.  (5)>  and 
b is  the  radius  of  the  conducting  cylinder.  Equation  (7)  indicates 
that  the  density  profile  has  sharp  radial  boundaries  at  ax  and  a2. 


The  azimuthal  current  density  is  given  by 


J§(r) 


qh  V°(r) 


o ^ r < ax  , 
ax  s r £ a2, 
a2  < r ^ b , 


where  Vg(r)  is  the  mean  azimuthal  velocity  of  the  p- layer,  i. e. , 

Vg(r)  - <vQ(r)>  = (l/n)Jv0/?d2v  = tur,  ax  * r s:  a2.  (9) 

Similarly,  it  can  be  shown  that 

V°  = <v  (r)>  = o. 

r r''  ' 

Thus,  the  p- layer  rotates  with  a constant  angular  velocity  id  around  its 
axis  of  symmetry  (rigid-rotor  equilibrium). 

The  vector  potential  corresponding  to  the  azimuthal  current  density 


of  Eq.  (6)  is  determined  from 


, o £ r < ax  , 


- (4tt/c)  qnruu,  ax  2 r ^ a2, 
o , ag  < r < b . 


The  solution  of  Eq.  (10)  subject  to  the  boundary  conditions  that  Ag  and 
B are  continuous  at  ax  and  a2  and  Jb2TTB  rdr  = TTb2Bo  (flux  conservation 


|k2r,  o 5 r 2 ax 

k2r  - TTnqu^r2  - a2)2/2cr,  ax  ^ r £ a2,  (11) 

k2r  - TTnqiu  [(a|  - a2)  (2I2  - a2  - af)]/2cr,  a2  < r ^ b, 
where  the  constant  k2  is  given  by 

k2  - B0/2  + (Trnqo)/2c)  (a2-af)  [2  - (af+al)/b2]  . (12) 

Using  Eqs.  (11)  and  (12)  and  the  relation  B°(r)  = ^ [rAg(r)], 

the  axial  magnetic  field  is  given  by 
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Bo  + 


. ld±di 


Sb2  J 


, ° ^ r 


ai 


B°(r)  = 


Bo  - * I,  [-^] 


, S2  < r ^ b , 


where  1^  = -a^  (af-af)  is  the  azimuthal  current  per  unit  length. 

Figure  2 shows  the  ratio  B°(r)/B0  as  a function  of  r/b. 

z 

The  inner  and  outer  radii  of  the  p- layer  are  determined  from  U(ai) 

= U(as)  = o,  where  the  function  U(r)  is  given  by  Eq.  (5)*  Substituting 
Eq.  (11)  into  Eq.  (5),  it  is  rather  straightforward  but  tedious  to 
show  that  the  radii  ax  and  a2  can  be  expressed  in  the  form 

l6ki  X^/fma^b2) 


a? 


and 


1 + {l  + l6(X*/b)2  ] ~(“ J + 1)  + 4kx ^ (mufV2)]  } 


(lb) 


a| 


b2  + 


SX*2  + 1)  + 2kx  / [mu^b2)] 


1 + {l  + l6(X*/b)2  [(^  + 1)  + 4kx / (m^b2)]  j172 


, (15) 


where  fJc  = qB o/mc,  u^i  a Mrnq 2/m  and  X*  = c/tu^.  Note  from  Eqs.  (l4) 
and  ( 15)  that 


. 1 + | 1 + i6(X*/b)2  ["(“J  + 1)  + 4k x / (mu.2b2)]  } , (16) 

and  since  a§  > a2,  the  denominator  in  these  equations  must  always  be 
positive.  Then,  in  order  for  ax  to  be  real,  kx  £ o. 

Radial  confinement  of  the  p- layer  requires  that  a§  ^ b2.  It  follows 
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then  from  ( 15)  that 


[(^  + 1)  - 2|kJ / (mu2b2)  ] £ o,  (17) 

where  |kil  “ -k1#  Furthermore,  since  (a2  - a2)  should  be  real,  the 
radical  in  Eq.  (16)  should  be  positive  or  zero,  i.e. , 

l6(A*/b)2  + l)  - 4|kJ / muPbz)  J ^ -1.  (18) 

Combining  Eqs.  (17)  and  (18),  it  is  obtained 


JOtl  ..  _j£_  * /nQ  + 1\sjal|4 
muFb2  16  X*2  y J wP b2 


(19) 


In  order  for  both  inequalities  to  be  satisfied  simultaneously. 


5!teA 


(20) 


irtJ^b2  ® 16\*2 

When  the  conducting  wall  is  absent,  i.e.,  b *♦  as,  equilibrium  exists, 
provided  u)  is  negative  for  positive  Q0  and  its  absolute  value  satisfied 
the  inequality 


o < |u)|  ^ Aq  . (21) 

Introducing  the  dimensionless  variable  \0  = ( 2 | kx | / m) 1 / 2 Qq1, 
Eq.  (19)  can  be  written  as 

/ ...  \ 

(22) 


^ (£> s 


where 


u>i  = [-!-(!  + 4X2/t>2)1^2]/2  , 


(23a) 


and 


u*3 


-l-fl  1-  f6Xp/b21  Cl  + (b/U»)g]P/g  ,„b, 

2[l  + (b/U*)*!  ' 1 V 1 


Figure  3 shows  % and  UJ3  as  a function  of  b/X0  for  several  values  of 
\*Ao.  For  small  values  of  \*A0  (high  ion  density)  there  is 
a small  range  of  b/\0,  in  which  equilibrium  exists  for  values 
of  I'Juj/Qo  » 1.  - In  general  such  equilibria 
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are  of  appreciable  thickness.  For  intermediate  values  of  X^Ao  (Fig 
the  frequency  of  rotation  is  restricted  to  o < |u)|  /C20  < 1.  In  the 
special  case  that  the  various  parameters  of  the  system  are  tuned  in 
such  a way  that  the  equality  of  Eq.  (20)  is  satisfied,  then  oh  °* 
and  the  frequency  of  rotation  is  given  by  Eq.  (23a).  Finally  (Fig.  3^) , 
for  large  values  of  X*A0  (small  ion  density),  the  normalized  frequency 
of  rotation  (co/Qo)  ^ id3  -1,  i.e.,  as  it  is  expected  all  the  ions 

rotate  with  a frequency  that  is  near  the  cyclotron  frequence  correspond- 
ing to  the  externally  applied  magnetic  field. 

It  is  easy  to  be  shown  from  Eq.  (16),  that  in  the  thin  layer  approx- 
imation the  frequency  of  rotation  (iu  is  given  by  an  expression 

similar  to  that  of  Eq.  (23b),  with  (X*)2  replaced  by  (X*)2/e,  where  e 
is  a positive  quantity  much  smaller  than  unity.  If  (b/4X*)  < 1,  the 
inequalities  of  Eq.  (22)  are  satisfied  and  the  thin  layer  frequency  of 
rotation  is  given  by 

(“thin700) 

By  contrast,  when  e ( b /4X*)  ^>>1  and  X0A*»1,  both  u&  and  (oj^^/Qq)  go 
to  zero,  i.e.  for  high  ion  density,  thin  layer  equilibria  exist  only 
with  very  low  rotational  frequency.  This  conclusion  is  similar  to  that 
of  LOVELACE  et  al.,  ( 1$ 74) . 

In  order  for  the  inequalities  of  Eq.  (22)  to  be  satisfied  (Jfe  s u^. 
Using  the  expressions  for  % and  ojq  from  Eq.  (23),  this  relation  is 
plotted  in  Fig.  4.  The  inequality  is  satisfied  only  in  the  upper  half 
of  the  figure.  For  fixed  values  of  XQ  and  b,  equilibria  exist  provided 
that  the  ion  density  in  the  layer  exceeds  a critical  value.  This 
clearly  demonstrates  the  importance  of  the  self  magnetic  field  in,  the 
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confinement  of  the  layer. 

III.  REVIEW  OF  INSTABILITIES 

The  stability  of  the  equilibrium  discussed  in  the  previous  section 
is  presently  under  investigation.  It  is  expected  that  results  will  be 
available  in  the  near  future.  In  the  mean  time,  in  order  to  obtain 
some  insight  about  the  stability  properties  of  the  rigid  rotor  equilib- 
ria, we  are  reviewing  several  unstable  modes  that  are  pertinent  to 
these  equilibria. 

Diocotron:  The  diocotron  mode  in  a space  charge  nonneutral  E- layer 

consisting  of  electrons  that  encircle  the  axis  of  symmetry  has  been 
studied  by  NOCENTINI  et  al.,  (I968).  They  found  that  the  sufficient 
condition  for  stability  is  the  absence  of  shear  in  the  velocity  of 
rotation.  Thus,  the  diocotron  mode  is  not  present  in  any  rigid  rotor 
equilibrium. 

Negative  mass:  The  stability  condition  for  a space  charge  neutral  E- 

layer  composed  of  relativistic  electrons  that  encircle  the  axis  of 


symmetry  and  confined  in  a magnetic  field  configuration  that  has  a 

r*  r 3B  -» 

weak  field  index  n I 1 = - ~ J is (LANDAU  et  al.,  1966) 

z 

^A/(2mccr0)J  > vgy/ct,  (stability),  (24) 
where  mo  is  the  electron  rest  mass,  g = c2L/(2nrQ),  *0  is  the  equilib- 


rium radius,  L is  the  inductance  of  the  beam  in  cgs  units,  A is  the 
momentum  spread  a “ l/(l-n)  - y“2,  v = Ne2/(2TTr0m0c2) , and  N is  the 
number  of  electrons  in  the  system. 

For  a strong,  nonrelativistic  p-layer,  both  the  field  index  and  the 
coefficient  a are  negative  and  the  stability  condition  [Eq.  (24)]  is 
always  satisfied,  i.e. , the  p-layer  is  stable. 
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However,  it  should  be  noticed  that  the  theoretical  work  of  Landau 
and  Neil  is  based  on  the  assumption  that  the  field  index  n is  small, 
which  is  not  probably  the  case  for  a field  reversing  p- layer.  Thus 
definite  conclusions  about  the  stability  of  a strong  p-layer  with 
respect  to  negative  mass  instability  cannot  be  drawn  before  recalculat- 
ing the  stability  condition  without  the  weak  field  index  assumption. 
Tearing  mode:  This  instability  is  not  of  any  great  concern,  because  it 

can  be  easily  stabilized  (MARX,  1968)  by  placing  a perfectly  conducting 
wall  sufficiently  close  to  the  charged  particle  layer. 

Precessional  mode:  It  has  been  predicted  theoretically  (FURTH,  1965) 

and  demonstrated  experimentally  ( CHRISTOFILOS  et  al.,  1 969)  in  the  case 
of  E- layers,  that  this  mode  is  present  only  in  weak  layers.  Further- 
more, it  has  been  shown  that  the  precessional  mode  can  be  stabilized 
by  a toroidal  magnetic  field  (CHRISTOFILOS  et  al.,  1972)  produced  by 
currents  flowing  along  an  axial  conductor  or  cantilever  or  by  a suffi- 
ciently strong  quadrupole  Ioffe  field  (WOODALL  et  al.,  1975)*  Thus,  it 
is  expected  this  mode  i .3  only  present  whenever  the  layer  strength  is 
below  a critical  level,  i.e.,  it  is  only  of  small,  if  any,  importance 
for  field  reversing  proton  layers. 
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IV.  CON CLUS IONS 


The  number  of  ions  required  to  form  a field  reversing  ring  or  layer 
is  approximately  by  the  factor  (tn./m  ) higher  than  the  number  of 
electrons  required  to  form  an  electron  ring  or  layer  of  the  same 
dimensions.  In  addition,  the  current  density  of  presently  available 
ion  sources  is  lower  than  those  of  electron  sources.  Therefore,  field 
reversing  ion  rings  or  layers  can  be  realized  by  using  hollow  ion  beams 
of  considerable  dimensions  and  passing  them  through  a magnetic  cusp. 
Since  the  canonical  angular  momentum  is  a function  of  radius,  such 
methods  of  injection  correspond  to  an  appreciable  spread  in  the  canon- 
ical angular  momentum  of  the  injected  ions.  Furthermore,  the  force 
that  an  ion  sees  as  it  passes  through  the  cusp  is  approximately  pro- 
portional to  its  radial  distance  and  thus  an  appreciable  spread  in  the 
energy  transverse  to  the  magnetic  field  lines  will  be  present.  There- 
fore, only  equilibria  that  are  based  on  distributions  functions  that 
include  spread  in  Pg  and  E^  are  pertinent  to  the  present  methods  of 
injections  (KAPETANAKOS  et  al.,  1976)  of  intense,  pulsed  ion  beams.  In 
this  respect,  our  selection  of  the  distribution  function  is  very  real- 
istic. However,  it  is  not  presently  clear  which  must  be  the  exact 
injection  conditions  in  order  to  form  an  ion  layer  in  which  all  the 
particles  have  the  same  H-u)Pg.  But  even  if  such  conditions  were  known, 
transient  effects  during  the  buildup  of  the  layer  could  drastically 
change  the  situation. 

Although  the  assumption  that  all  the  ions  have  the  same  H-xPg  is 
questionable,  it  has  the  important  advantage  that  it  results  to  a 


i 
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simple  distribution  function  that  allows  explicit  expressions  to  be 
obtained  for  the  particle  density,  current  density,  magnetic  field  and 
radii  of  the  layer.  Such  a distribution  function  is  very  convenient  for 
studying  the  stability  properties  of  the  layer. 
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